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[j. • Recent Monte Carlo simulation studies of a Lennard- Jones fluid confined to a mesoscopic slit- 

pore have reported evidence for "critical depletion" in the pore local number density near the 
liquid- vapour critical point, fn this note we demonstrate that the observed depletion effect is in 
fact a simulation artifact arising from small systematic errors associated with the use of long range 
corrections for the potential truncation. Owing to the large near-critical compressibility, these errors 

r^^ • lead to significant changes in the pore local number density. We suggest ways of avoiding similar 

^ ' problems in future studies of confined fluids. 
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I. INTRODUCTION 



-)— > 

In recent papers m-M , one of us has reported grand canonical Monte Carlo simulation studies of a Lennard- Jones 
I fluid confined between two structureless attractive walls arranged in a slit-pore geometry. The behaviour of the 

number density profile across the pore, p{z), was studied for various values of the thermodynamic parameters, namely 
the chemical potential, fj,, and temperature, T. At certain values of /i and T (apparently close to those of the bulk 
fj liquid-vapour critical point) it was observed that the average local density in the middle of the pore fell markedly 
below the value obtained in a fully periodic simulation performed at the same fj, and T. These findings were used 
to argue in favour of the existence of a generic "critical depletion" phenomenon, namely the proposed tendency of a 
critical fluid to be expelled by a conflning medium, even when the confining walls strongly attract the fluid particles 
#-/-\ . Such a scenario is supported by experimental findings for SFg adsorbed in mesoporous materials ||3|,y , for which 

^-H ' a dramatic reduction in adsorption was observed as the bulk critical temperature was approached from above along 
CsJ , the critical isochore. 

04 ■ In this note we point out that the apparent critical depletion reported in references 0-0] is actually a simulation 

^p [ artifact arising from systematic errors associated with the corrections applied to the configurational energy to compen- 
sate for the truncation of the interparticle potential. Using new simulations, we show that if one chooses a sufficiently 
large truncation distance or alternatively avoids the use of truncation corrections altogether, then the depletion effect 
jrt ■ disappears. 

^ ■ II. SIMULATION DETAILS AND RESULTS 

^ ' The simulation arrangement and procedure employed in this work are the same as those described in refs. [|l|,0, 

and accordingly we merely summarise the principal features. 

Grand canonical Monte Carlo simulations g were performed for a Lennard-Joncs fluid, having an interparticle 
potential of the form: 



UL.j{r) = 4e 
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where e and a are respectively the Lennard-Jones well depth and scale parameters. Two distinct geometries were 
studied: 

1. A fully periodic system. 

2. A slit-pore geometry, in which the fluid is confined between two parallel structureless walls, having periodic 
boundary conditions in the directions parallel to the walls. 



In the latter case, the walls were taken to exert a potential on the fluid particles of the form: 
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where / is a parameter that tunes the strength of the wall-fluid interactions relative to those of the fluid interparticle 
interactions. 

As in the previous studies of this system |3j3], the reduced temperature was set to the value T = 1.36, believed to 
be close to the bulk critical temperature. The chemical potential fj, of the periodic system was then tuned until the 
equilibrium density reached the value p = 0.365, believed to be close to the bulk critical density. The resulting value 
of /i was then fed into a simulation of the slit-pore system at the same temperature and with the choice / — 0.9836. 
In both the periodic and slit-pore arrangements, the Lennard- Jones interparticle potential was truncated at some 
radius and a compensating correction applied to the conflgurational energy. For the periodic system, this correction 
was calculated in the standard fashion by assuming a spherical cutoff surface of radius Vc centred on each particle, 
combined with a uniform density approximation (UDA) for r > r^. This yields for the energy correction per particle: 
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drr'^ULj{r) = -rTrpecr^ 
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where p — {N)/V is the average number density of the system. 

In the case of the slit-pore system, the truncation correction used was that given by Schoen et al. [p| , which assumes 
a cylindrical cutoff surface (on whose principal axis each particle lies) extending across the whole width D of the pore, 
i.e. such that the cylinder ends coincide with the pore walls. This yields 



^pore 



■ tan 



Sr 



(2.4) 



where Sc is the radius of the cutoff cylinder. 

We have studied the effect of the cylindrical cutoff radius Sc on the density proflle p{z) of the slit-pore system. For 
each choice of Sc, the p. value employed in the simulation was that yielding an average density p — 0.365 in a periodic 
system of the same dimensions and with Vc — Sc- For small cylindrical cutoffs (sc — 3.5ct), figure [^ shows that the 
local density in the pore middle is depleted with respect to the density of the periodic system (dashed line) at the 
same T,fj,. This is the result reported in refs. ||l|-^. However, new results for a larger choice of the cylindrical cutoff 
(also included in flgureH^) show that this depletion reduces as Sc is increased, and in fact vanishes for Sc -^ 5. Oct. This 
dependence of the depletion on the choice of Sc was missed in the previous studies [|l|,|[ . 

We have also performed simulations in which we dispense with the use of cutoff corrections altogether and simply 
simulate a system of particles interacting via a truncated Lennard- Jones potential. The results (figure 0), exhibit no 
sign of a density depletion in the pore middle with respect to the periodic system. 



III. DISCUSSION AND CONCLUSIONS 



The dependence of the density depletion on the choice of the cylindrical cutoff Sc (figure |]) points to a breakdown 
of the uniform density approximation invoked in the derivation of the truncation correction for the internal energy, 
eq. 2.4. This approximation assumes that the number density outwith the cutoff surface is uniform, having the average 
system density p = {N)/Y. However, figures |]andg show that for a slit system, p{z) exhibits considerable structure 
across the pore, especially close to the walls. Accordingly, one must expect some measure of systematic error to be 
associated with eq. ^^. Tests show that for the choice of cutoff Sc = 3.5(7 employed in references ||l|-H, this error is 



very small, so that in most circumstances eq. 2.4 represents a good approximation. 



It seems, however, that in the critical region, even a very small error in the truncation correction can lead to large 
effects on the local pore number density. The reason for this is the large near-critical compressibility, reffected in 
the fact that near Tc, isotherms of p{p) become very flat (see eg. figure 3 of H]). Since the error in the truncation 
correction acts rather like a shift in the bulk (chemical potential) field with respect to the periodic system, large 
alterations to the local pore density may result. This is in accord with the observation ||l|,^ that the depletion is large 
close to the critical point, but diminishes as one moves to higher temperatures along the critical isochore. 

To avoid similar problems arising in future studies of the effects of confinement on near critical ffuids, it would seem 
wise to adopt one of the following strategies: 



Employ a very large value for the truncation range and test for any systematic dependence of results on its 
value. This is clearly a very computationally intensive solution. 



2. Employ a truncated potential and dispense with corrections altogether, 
but causes complications if one wishes to model real substances. 



Such a system is clearly well defined 



3. Employ a cut and shifted potential which tends smoothly to zero at the cutoff. 

In summary, we have demonstrated that the apparent critical depletion reported in |l|-|^ was actually an artifact 
arising from systematic errors in the energy correction for the tail truncation in the slit-pore geometry. Although 
numerically small, these errors can (in the critical region) strongly influence the fluid local number density of the 
confined system compared to a periodic system at the same temperature and chemical potential. Thus our findings 
underline the care that must be taken when implementing any sort of truncation corrections for near critical fluid 
models. 

Further simulation studies of the near-critical properties of a confined fiuid are in progress and a detailed account 
of these and their implications for theory and experiment on critical depletion will be presented in a later paper. 
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FIG. 1. The number density profile p{z) (corresponding to a sht-pore of width D — 15a) at two values of the cylindrical 
cutoff radius. In each case, the reduced temperature is T = 1.36 and the chemical potential /j. used is that which yields an 
average density p — 0.365 (horizontal line) in a periodic system of linear size L = 15a^ at the same T. The figure shows that 
as the cylindrical cutoff radius is increased, the depletion disappears. 




FIG. 2. The number density profile p{z) (corresponding to a slit-pore of width D — 21a) for a truncated potential with no 
tail correction applied. In each case the reduced temperature is T = 1.36 and the chemical potential /i used is that which yields 
an average density p — 0.365 (horizontal line) in a periodic system of linear size L — 21(t^ at the same T. 



